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Abstract. The role of projectors associated with Poisson brackets of constrained 
Hamiltonian systems is analyzed. Projectors act in two instances in a bracket: in 
the explicit dependence on the variables and in the computation of the functional 
derivatives. The role of these projectors is investigated by using Dirac's theory of 
constrained Hamiltonian systems. Results are illustrated by three examples taken 
from plasma physics: magnetohydrodynamics, the Vlasov-Maxwell system, and the 
linear two-species Vlasov system with quasineutrality. 
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1. Introduction 

We consider an arbitrary Poisson bracket of a Poisson algebra of functionals of field 
variables given by 



where x G M", M" — )■ M"^, and Fy,-J-G-^ = J*-' G-^j with repeated indices summed. 
By Poisson algebra we mean a Lie algebra realization on functionals with an associative 
product of functionals that satisfies the Leibniz law. Also, we assume that the resulting 
equations of motion given by x = for some Hamiltonian functional 

possess a conservation law = 0, where Q is a functional of the field variables and 
their derivatives. Here we address the specific case where these conservation laws are 
obtained regardless of the choice of Hamiltonian if, so Q = is an intrinsic property of 
the bracket of the Poisson algebra. 

There are two ways to define such a constrained Poisson algebra. The usual way is 
to place a restriction on the set of field variables x the Poisson algebra. However, this 
definition raises the question of how to appropriately compute the constrained functional 
derivatives F^. The second way is to define a Poisson algebra that does not include any 
constraint on the field variables and, consequently, there is no ambiguity in defining 
the functional derivatives - conservation laws such as Q = take the form of Casimir 
invariants. 

In this article we investigate the links between these two ways of defining 
constrained Hamiltonian structures, and we propose a way to lift Poisson structures 
defined via the constrained field variables approach to ones that have the constraints 
as Casimir invariants. As can be expected, the difficulty resides in assuring the 
validity of the Jacobi identity. If we keep the same Poisson bracket but extended 
to the bigger algebra (the one without any constraint on the field variables), then in 
general, the Poisson structure is only obtained when the constraint is satisfied, i.e., 
the Jacobi identity is satisfied conditionally when = 0. It turns out that one can 
remedy this limitation by modifying the bracket with the inclusion of suitable projectors 
that leave the functional derivatives unconstrained and guarantee the Jacobi identity 
unconditionally. We identify such projectors acting on the functional derivatives and on 
the explicit dependence of the bracket on x- We discuss the various choices of projectors 
and highlight a particularly relevant one obtained from Dirac's theory of constrained 
Hamiltonian systems. 

In order to illustrate our purpose, consider the relatively simple and common 
example, the vorticity equation of a compressible or incompressible fiuid in R^. The 
vorticity a; = V x v, with v the velocity field, satisfies 




(1) 




(2) 



In terms of a commonly used Poisson bracket (see, e.g., Ref. [T]), 




(3) 
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Eq. ([2]) has the from F = {F,H}q with the Hamiltonian H = J (Pxv'^/2. However, if 
one forgets about the constraint on the vector fields u) or if one wants to lift the algebra 
of functionals of divergence-free uj to the algebra of functionals of any vector field u), 
then the bracket ([3]) does not satisfy the Jacobi identity. This is easily seen by the 
following counter example: 

Fi = - d^x (jj ■ ^y^ , F2 = - / d'^x u) ■ y z"^ , F3 = / d^x ■ zx , 



2.1 " ' " 2 

which yields, 



{Fi, {F2, FMo+ 0=- j d'xuj- V{yz) ^ . 



Evidently, the bracket (E]) satisfies the Jacobi identity only if V - a; = 0. We refer to such 
Poisson brackets that only satisfy the Jacobi identity conditionally as tainted brackets. 
One of the questions we address in this article is how to correct a tainted bracket so 
that it satisfies the Jacobi identity unconditionally. For this particular example, the 
correction is obtained by inserting a projection operator, following Ref. [2], given by 
V^ = l- VA-^V-, so that it defines a new bracket 

{F, G} = j d\ (Vuj) ■ (V X F^) X (V X G^). 

It is rather straightforward (see Ref. [2]) to show that this bracket satisfies the Jacobi 
identity unconditionally. We notice that V ■ is a Casimir invariant of the modified 
bracket, i.e. {V ■ a;, G} = for any functional G. 

As mentioned above, projectors are not only useful to lift algebras so as to 
satisfy the Jacobi identity, they are also involved in the way functional derivatives are 
computed when the field variables are constrained. As an illustration, we consider the 
incompressible Euler equation for the velocity field v(x, t). 



-V ■ Vv - VP. 



where P is determined by the constraint V ■ v = 0. This equation has a Hamiltonian 
structure [21111 El E] given by the Hamiltonian H[v] = f d^xv'^/2 and the Poisson bracket 



{F,G} = j rf^^xv- [Fv,Gv] 



where F^ are the functional derivatives of an observable F with respect to the field 
variable v and the Lie bracket [V, W]/, is given by 

[V, W]l = (W ■ V)V - (V • V)W. 

It should be noted that the incompressible Euler equation cannot be directly obtained 
from F = {F, H} using unconstrained functional derivatives F^ since V ■ v = 
would not be conserved by the flow. One way of correcting the bracket is to use an 
orthogonal projector P]. For divergence-free fields, this orthogonal projector is again 
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given by Pj_ = 1 — VA~^V- (see also Refs. [HE]). In other words, the constrained 
functional derivative must be computed such that it satisfies V ■ -Fv = 0. However, 
the fundamental reason for this constraint on the functional derivative is unclear, even 
though it yields the correct equation of motion. For a more general constraint =0, 
is it still the orthogonal projector that has to be used for the constrained functional 
derivatives? In addition, this projector is in general not unique. It therefore raises 
natural questions such as which is the most relevant projector and how is it obtained in 
a systematic way? 

In this article, we investigate two possible placements of a projectors: one is on 
the explicit dependence on the field variables, while the other is on the computation 
of the functional derivatives. We clarify the choice of the relevant projector by using 
Dirac's theory of constrained Hamiltonian systems. In order to prove the relevance of 
these projectors, we consider three examples taken from plasma physics. The first one is 
magnetohydrodynamics (MHD), both compressible and incompressible, the second one 
is the Vlasov- Maxwell system, and the third example involves semi-local constraints on 
linear Vlasov equations with two species. 

The goal of this paper is twofold: The first purpose is to show that using some 
relevant projector, the tainted brackets can be corrected such that the new brackets 
satisfy the Jacobi identity unconditionally. The second purpose is to connect these 
corrected brackets to the ones obtained from Dirac's theory of constrained Hamiltonian 
systems. 

2. Formulation of the general method 

2.1. Projected functional derivatives 

At the outset we assume that the bracket ([1]) is a Poisson bracket on the algebra of 
functionals of x, where x denotes a (i-tuple of fields such that Q[x](x) = and Q['x\ 
is function of x ^"^^ its derivatives. These fields will be referred to as Q-free fields. 
In this section, our aim is to get a corresponding Poisson bracket on the algebra of 
any functionals of Xi satisfying Q[x](x) = or not. The functional derivatives are 
defined in the following way: 



functional derivative F-^ from the unconstrained one F^ by the following equation: 




(4) 



for all Q-free 5%, which here means that Q ■ 6x = where Q is the Frechet derivative 
of Q defined by 





(5) 
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where now 5% is the constrained (Q-free) variation and 5x the unconstrained one. For 
the unconstrained variation 5%, we use a hnear operator V acting as Sx = V^^x such 
that QV^ = 0. Moreover, the range of this operator should be Ker Q and, in 
addition, should act as the identity on Ker Q. This is equivalent to requiring that 
V he a. projector. Consequently, this leads to a condition on the possible projectors V 
such that = VF^,, viz. 

Ker p = Rg Qt (6) 

Note that given this condition, Q[x](x) is a Casimir invariant that is naturally preserved 
by the flow. Still this projector is not unique. In the literature (see, e.g., Ref. [1]), the 
functional derivative is chosen such that QF-^ = 0, so that the projector satisfies QV = 0. 
This corresponds to the orthogonal projector 

V^ = 1-Q\QQ^)-'Q, (7) 

provided QQ^ is invertible on Rg Q. However it is not clear if it is the best choice for 
the projection. Other solutions satisfy 

Va_v = r^, 
vv± = V, 

which are needed in order to satisfy Eq. (j6]). Given a particular projector V the 
bracket ([I]) becomes 

{F, G}t = J d^x iVF^) ■ J(x) ■ (PG J, (8) 

where now the functional derivatives are the unconstrained ones. So we have released 
the constraint on the functional derivatives but, in general the Poisson bracket (|8]) does 
not satisfy the Jacobi identity for functionals of arbitrary Xi ones no longer restricted to 
Q-free fields, because JI(x) may give contributions that do not satisfy the Jacobi identify 
when Q[x\ 7^ 0. However, if the projector V does not depend on the field variables x, 
as is the case for the examples we deal with in this article, then a bracket that satisfies 
the Jacobi identity for all functionals of x, satisfying Q[x\ = or not, is given by 

{F, G} = J rx iVF^) ■ J(Px) ■ iVG^)- (9) 

In order to verify the Jacobi identity, we perform the change of variables Xp = T^X 
and Xq = X~ T^X so that bracket ([9]) formally becomes bracket ([H]) with xp instead 
of X- Since Xp is by definition Q-free, the Jacobi identity is satisfied. For the Poisson 
bracket (Q, we notice that Q[x](x) is a Casimir invariant, and that the equations of 
motion for xp ^'^'g identical to the ones given by the Poisson bracket ([I]) or ([8]). 



2.2. Dirac brackets 
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2.2.1. Local constraints In order to identify the most appropriate projector, we use 
Dirac's theory of constrained Hamiltonian systems [7]. We begin with the following 
good Poisson bracket: 

{F,G} = j rxF^-S{x)-G^ (10) 

and then impose the local constraint $(x) := Q[x](x) = 0, where as before Q[x](x) is a 
function of x(x) and its derivatives. The Dirac procedure begins with the computation 
of the matrix of Poisson brackets between the local constraints, 

C(x, x') = {<l>(x), <l>(x')} = QJQ^(5(x - x'). 

We set ^ := QJQt and we assume that this quantity is invertible. Then, the Dirac 
correction to the bracket (ITOl) is given by 




rf"xdV {F, $(x)}D(x, x'){$(x'), G}, 



where D(x, x') = A "'^(x(x))(5(x — x'). Since {F, $(x)} = — QJ- F-^, this contribution is 
equal to 

Therefore, the Dirac bracket is given by 

{F,G}, = J rxF^-Ux)-G^, (11) 

where 

J* = J- jqU-^qj. 

We notice that we only need to verify that A is invertible on the range of Q in order 
to define the Dirac bracket ( ITTi) . It is straightforward to verify that 1=,, is antisymmetric 
because A is antisymmetric. We notice that QJ* = (and therefore S*Q^ = 0). As 
a consequence, the constraint $ is a Casimir invariant. The Poisson brackets obtained 
by the Dirac procedure are Poisson brackets of the form ([8]) but untainted, i.e., they 
satisfy the Jacobi identity unconditionally even though they are not of the form ([9]) in 
general. This can be seen by considering a projector V as discussed in the previous 
section. Under the assumption that Ker V = Rg Q\ we deduce that J*(l — P) = 0, 
and consequently: 

With this equality, the Poisson bracket becomes 

{F, G}, = I (Tx (VF^) ■ Ux) ■ {VG^). 

The additional feature is that, a priori, the Poisson matrix is a function of both 
Vx and (1 — V)x. However, it is straightforward to check that (1 — V)x is a Casimir 
invariant. 
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Dirac's procedure shows that among the possible projectors V satisfying Eq. (E]), 
one turns out to be most convenient. The matrix J^, can be rewritten using the Dirac 
projector 

P. = 1-QU-^QJ, (12) 

as 

so that the Dirac bracket becomes the same as the original one (fTOj) with the exception 
that the functional derivatives are projected using the Dirac projector, 

{F, G}, = j d'^x iV^F^) ■ • (P.G^), (13) 

where we notice that the Poisson matrix is J and not J*. The main difference between 
the the orthogonal projector V±_ and the Dirac projector V* is that V±_ is a purely 
geometric object since it only depends on the constraints, and V* is a dynamical object 
since it involves the Poisson matrix. 

Remark: We observe that the matrix corresponding to the Dirac bracket has the 
following property: 

J* = vlw, = IV, = vli, 

i.e., the Dirac bracket can be rewritten from Eq. ( 1T3|) using only one Dirac projector, 
e.g., 

{F,G}, = J rxF^-J{x)-V.G^. 
As a result, the computation of the Dirac bracket is made easier. 

2.2.2. Semi-local constraints The calculation of Sec. 12.2.11 can be generalized to allow 
semi-local constraints in phase space. To this end we split the set of coordinates into 
two pieces, i.e., x = (xi,X2) where xi G M""™ and X2 G M™. The semi-local constraints 
are given by 

$(xi) = QW(xO = J rf™a:2QM(x), 

where Q[x](x) is a function of x(x) and its derivatives. The linear operator Q is defined 
by the linear operator associated with the function Q by 

Q = [ (rx2 Q. 



Since Q acting on a function of x is only a function of xi, the linear operator Q is 
defined by the equation 

r-'^xi Qx ■ w(xi) = I rxxi^) ■ q'w. 
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Consequently, Q is a linear operator acting on functions of xi as Q\ i.e., Q ly(xi) = 
Q}w{^i). In a manner similar to that of Sec. 12.2. ![ the computation of the Dirac bracket 
shows that the operator 

A = QIQ\ 

must be invertible. More explicitly, the linear operator A acts on functions of xi as 

^w(xi) = j d"'x2QIQ^wixi). 
The expression of the Dirac projector is given by 

n = 1 - Q^A-'QS, 

in a very similar way as the case of the local constraints. We notice that the linear 
operator A only needs to be invertible on Rg Q. Another important projector is the 
orthogonal projector given by 

As in the case of local constraints, these two projectors satisfy J* = V'^S^V, along with 
the two properties V±V^ = V± and V^V± = V^. In addition, the Dirac projector satisfies 

J* = vUv, = vU = JK- 

3. Example 1: magnetohydrodynamics 

3.1. Compressible magnetohydrodynamics 

A particularly interesting example is afforded by the Hamiltonian structure of 
magnetohydrodynamics. The equations for the velocity field v(x, t), the density p(x, t), 
the magnetic field B(x, t), and the entropy s(x, t) are given by 

P = -V-(pv), 

V = -V - Vv - p'^Vip^Up) + p~\V X B) X B, 
B = V X (v X B), 
s = —V ■ Vs, 

where U is the internal energy and Up here denotes the partial derivative of U with 
respect to p. The dynamical variables are p(x), v(x), B(x) and s(x) where x 
belongs to M^. The observables of the system are functionals of these vector fields, 
denoted generically by F{p, v, B, s). In these coordinates, this system has the following 
Hamiltonian 

d'xi-pv^ + pU{p,s) + — 
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There are two slightly different Poisson brackets that have been proposed in Refs. [SI El 
To] . A first one was given in Ref. [8], 

{F, G} = - j £x [FpV -G^ + F,- VGp - p-\V x v) ■ {F, x G^) 

+p-^Vs ■ (F.Gv - FvG,)] + {F, G}b, (14) 
where the magnetic part {F, G} s of the Poisson bracket is chosen as {F, G}b = {F, G} B,t 

{F, G}B,t = - j d^x p-^ (Fv ■ B X (V X Gb) - Gv ■ B X (V X F^)) . (15) 

It was pointed out in Ref. [lOj that this bracket satisfies the Jacobi identity only when 
V ■ B = 0, and also that V • B commutes with any other functionals, i.e., {F, V ■ B} = 
for all F (it is a Casimir-like property, even though we cannot call it a Casimir invariant 
since the Jacobi identity is only satisfied when V ■ B = 0). As was the case for the 
vorticity equation (l2l), the functional derivatives with respect to B must be divergence- 
free for coherence. However, we notice that here, since only V x Fb are involved in 
the expression of the magnetic part (fTT|) of the Poisson bracket, it does not make any 
difference whether Fb is divergence-free or not. 

In order to extend the definition of the Poisson bracket to functionals of any B, ones 
not necessarily divergence-free, a second Poisson bracket was proposed in Refs. P flO]. 
There the magnetic part of the Poisson bracket f|T^ was replaced by 

{F, G}b,i = - Jd^x [(p-^Fv • [VGb] - p-'G, ■ [VFb]) ■ B 

+B ■ ([V (p-'Fv)] ■ Gb - [V [p-'G^]] ■ Fb)] . 

Here the notation a ■ [M] ■ b is a scalar explicitly given by Ylij O'i^ij^j for ciny vectors a 
and b and any matrix [M]. It was shown that this bracket satisfies the Jacobi identity 
for all functionals of (p, v, s, B) regardless of the condition V ■ B = 0. The magnetic 
part of this Poisson bracket is rewritten as 

{F, G]b,i = - j d^xp-^ [Fv • B X (V X Gb) - Gv ■ B X (V X Fb)] 

+ j d^x p-i V ■ B (Fv ■ Gb - Fb ■ G^) . (16) 

The first line of the above bracket corresponds to the Poisson bracket introduced in 
Ref. [Hj [see Eq. (fTSl) ]. With the additional terms (proportional to V ■ B) the Jacobi 
identity is unconditionally satisfied for any functionals of (p, v, s,B). However, a 
property of the bracket f|T^ with the magnetic part f|T5|) has been lost, V ■ B does 
not Poisson-commute with any functional, so it is not a Casimir invariant. 

In order to have both the Jacobi identity unconditionally satisfied and V ■ B a 
Casimir invariant, we apply the prescription ([9]) on the magnetic part (|T5l) . At every 
instance in the Poisson bracket where B is explicitly mentioned, we replace B with 
B = B — VA~^V ■ B. The magnetic part becomes 

{F,G}b = - j d'^xp-' (Fv- (B X (V X Gb) -Gv- (B x (V x Fb))) , 
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and it is rewritten as 

{F,G}b = - j d^xp'^ (Fv • (B X (V X Gb) - Gv ■ (B x (V x Fb))) 

+ j d^xV -B A- V • {p-^F^ X (V X Gb) - p"^Gv x (V x Fb)) . (17) 

Here we notice that the correction term still contains terms proportional to V ■ B but 
is different from the one in Eq. ( fT6l) . The main difference is that V ■ B is not a Casimir 
invariant for the Poisson bracket ( IT6l) whereas it is one for the Poisson bracket ( IT7|) since 
it only involves terms like V x Gb- 



3.2. Incompressible magnetohydrodynamics 



For incompressible MHD we begin with the equations for compressible magnetohydrody- 
namics from Sec. I3.1l and apply constraints. The Poisson bracket given by Eqs. f|T4|) -f|ri 
is of the form fITOl) with 



/ 



V 








-V- 

— p"^(V X v)x 
-V X (p^^Bx) 
-p-^Vs 







-p-^B X (Vx; 








\ 

p-^Vs 





We impose the following local constraints on the field variables % = (p, v, B, s), 

Q[x](x) = (p,V-v). 

The reduction to incompressible MHD using Dirac's theory has already been done in 
Ref . PI and the reduction to incompressible Euler equation in Refs. [HI [12] . Here we 
propose a more compact way to present this reduction using the operators introduced 
in the previous sections. The expressions of the intermediate operators are 



Q 



A 








A 



A 



-1 



-A V- (p-i(V X v) X V) 

A-^V- (p-i(V X v) X V) 
A-i 



-A-i 




The orthogonal projector is given by Eq. ([7]) and its expression is 



P^F^ = (0,Fv,i^B,i^s), 
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where = — VA ■ F^. The Dirac projector, computed from the Poisson 
bracket f lT^ where B has been replaced by B = B — VA^^V ■ B, is given by 

T'^ifF-^ i.-F*1 F'Vl -^B; Ps)i 

where 

= A'V ■ {p-^ ((V X v) X - B X (V X Fb) - F,Vs)) . 

We notice that the two projectors differ in the first component. Even though the 
two projectors V±_ and are different, both of these projectors satisfy the equation 
J* = P^JP, which is always the case for the Dirac projector but not true in general 
for the orthogonal projector. Actually any projector VF^^ = (F*(Fv, Fb, Fg), FV, Fb, Fg) 
satisfies J^. = V'^SV for any function F^. The first component is thus irrelevant, and 
consequently the orthogonal projector is the simplest projector to be used for constrained 
functional derivatives. From this projector, we compute the Dirac bracket from Eq. ( !T3|) . 
and it gives the same bracket as that produced in Ref. [2]: 

{F, G}, = J £xp-' ((V X v) ■ (FV X Gv) - Vs ■ (F.Gv - F^Gs) 
+B ■ (Fv X (V X Gb) + (V X Fb) x G^)) , 

where F^ = F^ - VA^V ■ F^. 



4. Example 2: Vlasov-Maxwell equations 

4.1. Vlasov-Maxwell modified bracket as a Dirac bracket 

As a second example, we consider the Vlasov-Maxwell equations for the distribution 
of charged particles in phase space /(x, v,t) and the electromagnetic fields E(x, t) and 
B(x, t) given by 

/ = _v.V/-(E + vxB)-9v/, 
E = V X B- J, 
B = -V X E, 

where J = / d^vvf. The Hamiltonian of this system is given by 

H = Jd^zf^ + ld^x^^±^, 

where we denote z = (x, v). The Poisson bracket between two functionals of /(x, v), 
E(x) and B(x) is given by 

{F,G}t= j d'zf{[Ff,Gfl + [Ff,Gf]B + GE-d^Ff-FE-d^Gf) 

+ y c/=^x (Fe ■ V X Gb - V X Fb ■ Ge) , (18) 
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where the two brackets [■, -jc and [■, -Ib are defined by 

[f,g], = Vf -d^g-d^f -Vg, (19) 
[f,g]^ = B-{d^fxd^g). (20) 

The Poisson bracket (|T8|l was given in Ref. [13] based on an earher work of Ref. [14J (see 
also Ref. [I5j). It was pointed out in Ref. [lO] that the Poisson bracket ( |T8i) only satisfies 
the Jacobi identity when V • B = 0, which is to say that it does not satisfy the Jacobi 
identity for arbitrary functionals of (/, E, B). This problem is actually already present in 
the Lagrangian description (for the dynamics of charged particles) since [■, -Ic + l', ■]b only 
satisfies the Jacobi identity for functions B such that V ■ B = 0, whereas, individually, 
[■, -Jc and [•, -Ib satisfy the Jacobi identity for an arbitrary function B. 

In order to remedy this problem, we modify the bracket [■, -Jb to take the form of 

o, 

[/,^?]b, = (B - VA-V ■ B) • {d^f X d^g). 

With this modified gyrobracket, we readily check that [■, -jc + [■, ■]bp satisfies the Jacobi 
identity. Next, we consider the modified Poisson bracket (|T8|) obtained by replacing 
[■, -Jb by [■, -Jbp, i.e., we consider the Poisson bracket 

{F,G}vM= j d'zf{[Ff,Gfl + [Ff,Gf]B,+GE-d^Ff-FE-d^Gf) 

+ j d'x (Fe-V xGb-V X Fb- Ge) , (21) 

which satisfies the Jacobi identity unconditionally. This follows from the change of 
variable Bp = B - VA'^V ■ B and Bn = VA^^V ■ B where it should be noted that 



V X Gb = V X Gb 



p ' 



since the operator V = 1 — VA^^V- satisfies Wx = Vx. 

Here it should be noticed that V ■ B is a Casimir invariant for the Poisson 
bracket (l2Ti) . The untainted form of the Vlasov-Maxwell bracket ( 1211) gives the 
Hamiltonian structure of the Vlasov-Maxwell equations in terms of physical fields 
without introducing the vector potential, i.e., without the restriction of V ■ B = 0. 

In order to realize the link between brackets defined using projectors and Dirac 
brackets, we show below that the Poisson bracket fl2T|) is a Dirac bracket of some parent 
bracket obtained using two constraints which, by definition, are Casimir invariants of 
the bracket (|2T|) 

Q[/,E,B](x) = (V-E-p,V-B), 

where p = J d^v f. As expected there is an infinite number of solutions for the parent 
bracket. A family of solutions is given by 



{F, G} = {F, G}vM + J d'x (V ■ FbPV ■ Ge - V ■ F^V^V ■ Gb) 



(22) 



where P is a linear operator independent of the field variables, so that the Jacobi identity 
is guaranteed by Morrison's lemma of Ref. . This statement uses the fact that the 
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Vlasov-Maxwell bracket has been made untainted; it would not be true if the original 
tainted Vlasov-Maxwell bracket (|T8l) was considered instead of the Poisson bracket f l2T]) . 

Now, if we apply the Dirac procedure on the extended Poisson bracket f l22l) with 
the primary constraint V • E — p, we get the secondary constraint V ■ B, and the reduced 
Dirac bracket is obtained from J.,, = P^JP* where is the Dirac projector (fT2l) . The 
Dirac projector can be explicitly computed. However, in order to further simplify the 
computation of the Dirac bracket, we use the orthogonal projector since, as in the case of 
incompressible MHD (see Sec. 13. 2p . it satisfies the same relation as the Dirac projector, 
i.e., J, = VUV, = V[IV±, where 

V±F^ = (Fb - VA-V ■ Fb, Fe, Ff). 

This implies the expected result that the Vlasov-Maxwell bracket (l2Ti) is the Dirac 
bracket of the bracket (|22|) with Dirac constraints (V ■ E — p, V ■ B). 

With the extended bracket ( 122|) . the Casimirs (V-E— p, V-B) of the Vlasov-Maxwell 
system now have dynamics given by 

^(V-E-p) = AI)tv.B, 

d 

_V-B = -APV-E. 

dt 

These equations suggest two particularly interesting choices for our still undetermined 
operator V. Defining D = A^^ gives to V ■ E — p and V ■ B the dynamics of stationary 
waves when p = 0, whereas defining T) = (— A)^^/^ gives them the dynamics of 
propagating waves. We note that these operators always act on divergences of vector 
fields. 

Remark: As a side note, we point out that the choice of P = (— A)~^/^ 
naturally exhibits the operator V* := V(— A)^^/^V- which corresponds to Vx for the 
compressible part of a vector field. Indeed, the operator 

-V* A-V* = VA-V-, 

is the orthogonal projector onto the kernel of Vx, just as —V x A^^Vx = 1 — VA^^V- 
is the complementary projector onto the kernel of V*. With this choice, the resulting 
dynamical equations associated with the Poisson bracket ( 122|) for the solenoidal and the 
compressible parts of the electromagnetic fields become independent and similar: 

UEs = —Js, DEc = —Jc, 

DBs = V X J5, DBc = V * Jc, 

where □ is the d'Alembert operator □ = d'^/dt'^ — A and il^s is the solenoidal part of the 
vector field ij), i.e., ips = —V x A~^V x -0 = (1 — VA~^V) -ip and ipc is its compressible 
part, which is i/'c = —V * A^^V = VA^^V • if). In the absence of matter, the fields 
'05 and V'c propagate as independent free waves. 
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In order to obtain Vlasov-Poisson equations from the Vlasov-Maxwell equations we 
impose two constraints: 

Q[/,E,B](x) = (B-Bo(x),VxE), 

where Bq is a non-uniform background magnetic field. The operators Q and are 
given by 

^ V- 

1 



Q 



and 



/ 0\ 



\ 



V- 

1 



The orthogonal projector V± given by Eq. is given by 



1 \ 
VA-^V- 
y 



Contrary to the orthogonal projector, the expression of the Dirac projector depends on 
the dynamics, and in particular on the Poisson matrix J which is given by 

/-[/,-] -d^f \ 
-/9v Vx 
-Vx 



where the small bracket [■, ■] is given by [■, •] = [■, -Jc + [■, -Jsp with these two brackets 
given by Eqs. ( fT9l) - (l20l) . The operator A is given by 



A 



(Vx)2 
-(Vx)2 



The operator A is not invertible, but one can give an expression for ^ by 
A-'Q = 



A-i(l - VA-^V-) 
-A-^Vx 



As a result, the Dirac projector is computed, 

/ 1 



1 

VA-^V- 

\^ -A-^V X /9v vA-iy- 
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We notice that both projectors V±_ and satisfy the equation 
Poisson matrix of the Vlasov-Poisson equations is given by 

/ 



V^V and the 



v 







-VA-iy ■ 9v/ \ 








It leads to the expression of the Poisson bracket, 

= j d'z f[Ff - A-'V ■ FE,Gf - A-'V ■ Ge]. 

Like for the incompressible MHD equations, even if the Dirac and orthogonal projectors 
are different, both of them can be used to compute the Dirac bracket from the Poisson 
matrix J, the orthogonal projector being slightly simpler and more straightforward to 
compute. 



5. Example 3: Quasi-neutrality as semi-local constraints 

In this section, we give an example of a set of physically relevant constraints where 
the orthogonal projector does not exist, and where the Dirac projector is a natural 
replacement for computing the constrained functional derivatives. We consider the 
following Vlasov equation with two species, ions and electrons, linearized about spatially 
homogeneous distribution functions aj(v) and ae(v). The equations for the phase space 
density fluctuation of ions /i(x, v) and electrons /e(x, v), are given by 

fs = -V ■ Vfs + Cs d^as ■ V0, 

where A0 = — J2s J d^v fs and = ±1 is the charge of the particles of each species 
s = i, e. 

The field variables are x(z) = (/i(z), /e(z)), which are functions of z = (x, v). The 
Poisson bracket, in this case, is defined by the Poisson matrix J given by 

[«i(v),-] 
[«o(v),- 

with [^^(v),^] = —dvOs ■ VG. According to Ref. Pj (see also Ref. [IS]), the 
corresponding Hamiltonian can be found as the quadratic functional corresponding 
to the second derivative of the Hamiltonian (plus Casimirs) of the nonlinear Vlasov- 
Poisson system, evaluated at fs = as(v). This is true for ion and electron densities 
close to equilibria that are isotropic in velocity (like a Maxwellian for instance). We 
impose the set of two semi-local constraints 

Q[X](X)= dh{f,-f,)J d3^vV(/i-/e)y 

The first component of the constraint is the quasi-neutrality. The second component is 
a secondary constraint associated with quasi-neutrality according to Dirac's theory of 
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constrained Hamiltonian systems [T7]. Since the constraints are linear with respect to 
the field variables, the operators Q and Q are given by 

The operator Q acts on functions of x only. In order to compute the Dirac projector, 
the matrix A = QSQ needs to be computed 



A 



(oi + ae)A 

(ai + ae)A 2(A + /3c)-VA 



with as = J d^vUs and f3s = J d^v\as- This operator is invertible and its inverse is 

A = _ _ fti + ftc 

Note that the operators A and A^^ act on functions of x and return a function of x. 
The Dirac projector V^. has the form 

A-^V- /-^g , ^ / [a;,.] -[a,,-] 

= 1 - ^ — d V {v + V - 2v) \ \ \ / ' 

ai + ae J y -[fti, -J [fte, -J 

where v = (A + /3e)/(ai + Oe)- Concerning the orthogonal projector, we note that QQ 
given by 



2 / ' 



does not exist since it is unbounded when it acts on functions of x. As a consequence, the 
orthogonal projector cannot be a solution for the computation of constrained functional 
derivatives. Here a convenient choice is afforded by the Dirac projector. 
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